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ABSTRACT 

Several aspects of the 4d imprint of the 5d bulk Weyl radiation are investi- 
gated within a recently proposed model solution. It is shown that the solution 
has a number of physically interesting properties. The constraints on the model 
imposed by combined measurements of rotation curve and lensing are discussed. 
A brief comparison with a well known scalar field model is also given. 

1 INTRODUCTION 
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Early observations led to the hypothesis that there could be large amounts 
of nonluminous matter hidden in the galactic haloes (Oort 1933; Zwicky 1933, 
1937). Later observations of flat rotation curves of spiral galaxies confirmed 
the hypothesis (Freeman 1970; Roberts & Rots 1973; Ostriker, Peebles & Yahill 
1974; Einasto, Kaasik & Saar 1974; Rubin, Thonnard & Ford 1978; Rubin, 
Roberts & Ford 1979; Sofue & Rubin 2001; Kochanek et al. 2005). Doppler 
emissions from stable circular orbits of neutral hydrogen clouds in the halo allow 
the measurement of tangential velocity Vtg{r) of the clouds treated as probe 
particles. According to Newton's laws, centrifugal acceleration v^g/r should 
balance the gravitational attraction GM(r)/r'^, which immediately gives = 
GM{r)/r. That is, one would expect a fall-off of Vfg{r) with r. Observations 
indicate that this is not the case: Vtg approximately levels off with r in the halo 
region. The only way to reconcile this result of observation is to hypothesize that 
the mass M(r) increases linearly with distance r. Luminous mass distribution in 
the galaxy does not follow this behavior. Hence the conclusion that there must 
be huge amounts of nonluminous matter hidden in the halo. This imsecn matter 
is given a technical name "dark matter" . Gravitational lensing measurements 
have further confirmed the presence of dark matter (Maoz 1994; Barnes et 
al. 1999; Cheng & Krauss 1999; Sofue & Rubin 2001; Trott & Webster 2002; 
Weinberg & Kamionkowski 2002; Kochanek & Schechter 2004; Smith et al. 
2005; Faber & Visser 2006; Metcaff & Silk 2007). Current estimates suggest 
that about 23% of matter in the whole universe consists of dark matter residing 
in the galactic haloes. 

Although the exact nature of dark matter is as yet unknown, several candi- 
dates for it have been proposed in the literature. One of the favored candidates 
is the Standard Cold Dark Matter of the so called SCDM paradigm (Efstathiou, 
Sutherland & Madox 1990; Pope et al. 2004). Despite its initial success, the 
current consensus seems to converge on the Lambda CDM or ACDM model that 
is related to the accelerating expansion of the universe (Tegmark et al. 2004 
a,b). 

Analytic halo models include the framework provided by scalar-tensor theo- 
ries. Scalar fields are important because they are predicted by supersymmetric 
unification theories. (See for instance, Ellis et al. 1998.) In particular, a proto- 
type of scalar-tensor theories, namely, the Brans-Dicke theory with scalar field 
(j) has the potentiality to explain a wide range of effects, from those in solar sys- 
tem (Weinberg, 1972; Bhadra, Sarkar & Nandi 2007) and gravitational lensing 
(Bhadra 2003; Nandi, Zhang & Zakharov, 2006; Sarkar & Bhadra 2006) to those 
arising from objects as exotic as wormholes (Agnese & La Camera 1995; Nandi, 
Evans & Islam 1997; Nandi et al. 1998; Bhadra & Sarkar 2005). Fay (2004) 
considered a larger class of scalar-tensor theories with a potential V{(f)), and 
with coupling parameter to = lv{4>), for the investigation of galactic dynamics. 
Many other halo models including suitable variants of scalar field theories also 
exist in the literature. See for instance, Bekenstein & Milgrom (1984), Sanders 
(1984, 1986), Soleng (1995), Matos, Guzman & Nufiez (2000), Peebles (2000), 
Matos & Guzman (2001), Nucamendi, Salgado & Sudarsky (2001); Mielke & 
Schunck (2002); Cabral-Rosetti et al. (2002), Lidsey, Matos & Urena-Lopez 
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(2002), Bharadwaj & Kar (2003); Arbcy, Lesgourgucs & Salati (2003). Sub- 
stantial halo models are provided by brane world theory (Mak & Harko 2004; 
Rahaman et al. 2008). 

The motivation for the branc world model, which we are considering in this 
article, comes from an entirely different but important direction. It comes from 
the consideration of higher dimensional spacetime because it can be argued 
that Ad is not big enough to fully accommodate the self interaction dynamics 
of gravity. In a recent article, Dadhich (2009) has given persuasive arguments 
based on flat space imbedding, self interaction and charge neutrality. But why 
only 5(i and not higher? The reason is that the conformally flat FRW universe 
is imbeddable only in 5d flat spacetime, which means that free gravity does 
not propagate farther than 4d because Weyl curvature vanishes for this FRW 
universe. The Randall-Simdrum (1999a, b) model provides the following setting: 
The 3— brane (that is, our 4d spacetime which confines matter and other gauge 
fields) has zero Weyl curvature and it bounds the 5d constant curvature bulk 
spacetime. To probe into fifth dimension, one can not rely on light signals as 
they do not propagate there. Therefore a possible probe has to rely only on 
propagation of gravity off the brane in an unfamiliar non-free manner (Dadhich 
2009). For this, one needs to devise a purely gravity experiment and a method 
to fathom gravity propagation in higher dimension. This being a formidable 
task, the usual approach is to look for measurable imprints coming from 5d into 
known gravitational configurations. The galactic halo, where gravity is the king, 
provides a natural large scale arena for this observation. It is quite probable 
that the effect of unseen dark matter is due to such imprints which contribute 
to the source term of Einstein's Ad field equations. Solving Einstein's equations, 
one then tries to work out the details of higher dimensional non-local effects. 

Static spherically symmetric exterior solutions of the brane world model 
have been derived in the literature (Dadhich, Maartens, Papadopoulos & Reza- 
nia 2000; Germani & Maartens 2001; Casadio, Fabbri & Mazzacurati 2002; 
Visser & Wiltshire 2003; Harko & Mak 2004; Mak & Harko 2004; Creek et al. 
2006; Viznyuk & Shtanov 2007). While all of these works are motivated by one 
or the other relevant considerations, the solution by Rahaman et al. (2008) is 
based on the simplest supposition of rotation curves that connects dark matter 
with the brane world model in a straightforward way. Whatever be the analytic 
model, there must be a way to contrast its predictions with actual measure- 
ments. The key point is that one does not directly measure the metric functions 
but indirectly measures gravitational potentials and masses from rotation curve 
and lensing observations. Bharadwaj & Kar (2003) advocated that a combina- 
tion of rotation curves and lensing measurements could be used to determine 
the equation of state of the halo fluid, although their model was restricted to 
specific form of flat rotation curve as well as guesses on equations of state. Faber 
and Visser (2006) have shown how, in the first post Newtonian approximation, 
the combined measurements of rotation curves and gravitational lensing allow 
inferences about the mass and pressure profile of the galactic halo as well as its 
equation of state. 

In this paper, we wish to explore various features of the brane world solution 
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by Rahaman ct al. (2008) focusing in particular on the constraints imposed 
by combined measurements. Our general conclusion is that the model has a 
number of interesting properties. Nevertheless, it is too premature to say if the 
model is obscrvationally supported. The relevant measurement constraints on 
the model are worked out. The paper is organized as follows: In this section, we 
have already delineated the motivation for the brane world model. In Sec. 2, we 
briefly outline the fleld equations and the solution under consideration for later 
use. Sec. 3 discusses the nature of the galactic fluid while Sec. 4 shows stability 
of circular orbits. In Sec. 5, we establish the attractive nature of gravity in the 
halo arguing from two different viewpoints and in Sec. 6 we discuss observational 
constraints. In Sec. 7 we make a brief comparison with the central features of a 
scalar field model. Sec. 8 contains the conclusions. 

2 FIELD EQUATIONS AND THE SOLUTION 

The 5d bulk field equations with negative vacuum energy A5 can be reduced 
to "effective" 4d gravitational field equations on the 3d brane. The governing 
equations are due to Shiromizu, Maeda & Sasaki (2000): 

Gij = Rij — —QijR = —Agij + k\Tij + k^Sij — Eij, (1) 

where the Roman indices nm from to 3, fc^ = 87rG5 is the 5d gravitational 
constant, A is the 4d cosmological constant given by A = fc|(A5 + fc|A^/6), the 
4d coupling constant k4 is fcl = and Xb is the vacuum energy on the 

brane. The stress E^j is the non-local effect from the bulk given by (Dadhich et 
al. 2000) 



where k = kc,/k4, hij = gij +UiUj is the projection tensor, U = k'^EijU^u^ is the 
so-called "dark radiation" energy density, Qi is a spatial vector while is a 
spatial, trace free symmetric tensor. The E^j satisfies, in virtue of the Bianchi 
identities, the conservation law 

EU = klS^^i (3) 

where semicolon denotes covariant derivative with respect to gij. The simplest 
case is the case of vacuum (absence of ordinary matter) so that = =^ Sij = 
which, in the paradigm of the study here, implies that we are interested in the 
region away from the galactic core, where the dark matter strongly dominates 
over the ordinary matter. In static vacuum, we have Qi = and in the comoving 
orthonormal frame, we have = 6q, hij =diag(0, 1, 1, 1). In this case, the 
constraint Eq (3) suggests a solution (Dadhich et al. 2000) 

P^j = P[nr, - (4) 
where r, is a unit radial vector and P is the so-called "dark pressure" . 
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Let us consider a spherically symmetric line element 



(5) 



where v and A are the metric potentials and are function of the space coordinate 
r only, such that v = u{r) and A = A(r). For this symmetry, Pij = P{r)[rirj — 
^hij]. In static vacuum {Tij = 0, Qi = 0) with the trace free source term El = 
-R- = 0, and with the unit radial vector = (0, 1, 0, 0), the field equations 
(1) yield the following (Mak & Harko 2004): 
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where U = U{r), P = P(r) and primes denote differentiation with respect to r. 
[Eq.(7) of Rahaman et al. (2008) has been properly rearranged in Eq.(8) above, 
which reveals the correct Ricci component on the left hand side] 

Assuming the known flat rotation curve condition v^g 
(Chandrasekhar 1983), Rahaman et al. (2008) obtained a solution as follows: 



^-^^^ = constant 



(10) 



A{r) 



2 D 
(2 + i)a ^ ^ 



(11) 



where I = 2{vtg) , Bq > and D are arbitrary constants and a = (2 + Z + 
y)/(2+^). They showed that the mass function M(r) = J ?7r^(ir increases 
linearly with r, which agrees with observation (Begeman 1989). However, this 
expression for mass is Newtonian and is applicable only so long as one can 
neglect pressure contributions although as yet there is no observational ground 
for such a neglect. To see the expressions for various types of mass distributions 
in terms of observational parameters, we wait till Sec. 6. Meanwhile, there are 
several other facets of the solution that need attention, which are addressed 
below. 

3 NATURE OF DARK RADIATION 

To throw some light on this issue, we identify the right hand side of the 
equations (6)- (8) with the standard spherically symmetric fluid distribution (p. 
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Pr, Pt, Pt) via Einstein's equations Gjj = —SttEjj in the observer's orthonormal 
rest frame ("). Using Eqs.(lO) and (11) in (6)-(8), we get: 



l{l + 2)r-^ Dr-°'l{l + l) 



+ 



1 + 4 



U + 2P 



Dr-°'{1'^ + 3/2 + 6/ + 4) - (/ - 2)r- 



/3 



(12) 



(13) 



— Gfn(n 



U~P 



'00 



a 



vv 



P{1 + 4)r-2 - 2L'r-"(/3 + 3P + 61 + 4) 



+ 12/ + 16 



/2 + 4Z + 12 

1 + 4 ' 



= 1^ + 21 + 4 



(14) 
(15) 



where we have chosen units on the right hand side of Eqs.(6)-(8) such that 
= 1 and redefined the fluid stresses (p, pr, pt, Pt) in terms of U and P. 
Eliminating U from Eqs.(13), (14), one can easily find the expression for P{r) 
as given in Rahaman et al. (2008) after their arbitrary constant of integration 
is set to zero. There is actually no room for a nonzero constant because once 
the metric functions are known, U and P follow from Eqs. (6)-(8). 

The interesting results we find are the following: While U is the same as p, 
the dark pressure P{r) is equal to neither p,, nor pt but P{r) = Pr{r) —pt{r) ^ 0. 
This indicates that the dark Weyl radiation imprint can not be characterized 
in the rest frame by a perfect fluid which requires Pr{r) = Pt{r)- So the perfect 
fluid version of dark radiation is that U{r) ^ but P(r) = 0. Looking at Eq.(2), 
we see that the imprint radiation in Ad vacuum is a dust-like fluid where Eij = 
—k^U{r) {uiUj — \hij ) . This is not the case with the solution under consideration 
since Pr{r) ^ Pt{r). This pressure anisotropy is a good feature of the solution 
from the point of view of exterior matching. Note that the solution cannot be 
matched to the Schwarzschild exterior metric at the boundary of the halo if the 
pressures were isotropic (Bharadwaj & Kar 2003). The Schwarzschild solution, 
although not the unique spherically symmetric vacuum when considering the 
brane-world model, is desirable for mimicking the far region (vacuum) of true 
Ad gravity, as supported by observation. We further see that the dark radiation 
fluid is not of exotic nature. The minimal condition for exoticity is that the Null 
Energy Condition (NEC) should be violated (Hochberg feVisser 1998). We find 
that the NEC is satisfied at all radii because p + Pr > (Fig.l). Transverse 
pressures are not included as they refer only to normal matter (Visser, Kar & 
Dadhich 2003). Even if we include them, we find p + Pr + 2pt > for all r. 

In the rest of the article, we assume the following numerical values for con- 
stants. Observations of the frequency shifts in the HI radiation show that, in 
the halo region, Vtg/c is nearly constant at a value 7 x 10^^ (Binney & Tremaine 
1987; Persic, Salucci & Stel 1996; Boriello & Salucci 2001). Thus I ~ lO'^. We 
also note that the solution is not used for small r, the inner core of the galaxy. 
Accordingly, to be consistent with observational facts, in all the figures in this 
paper we shall take Bq = 1, 1 = 0.000001, D = 0.00001 and fairly large distances 
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in Kpcs. 

4 STABILITY OF CIRCULAR ORBITS 

Defining the four velocity U" = of a test particle moving solely in 
the subspacc of the branc (and restricting ourselves to = 7r/2), the equation 



-mi can be cast in a Newtonian form 



(16) 



which gives 
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where the constants E and L, respectively, are the conserved relativistic energy 
and angular momentum per unit rest mass of the test particle. Circular orbits 
are defined by r = i? =constant so that ^ = and, additionally, ^ \r=R= 0. 
From these two conditions follow the conserved parameters: 
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and using it in V{R) = —E"^, we get 



E = ± 



2 - 



(20) 



The orbits will be stable if ^ 



r=R< and unstable if |r=j?.> 0. Putting 



the expressions for L and E in 
culations, the final result, viz.. 



\r=R, we obtain, after straightforward cal- 
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(21) 



Thus \r=R< so that circular orbits are stable in the model under consid- 
eration when D > 0. The latter condition is always satisfied since, from Eq. 
(11), we see that D has a dimension proportional to a power of radius R. The 
radius is always positive and so is D. 

5 ATTRACTION IN DARK RADIATION 

Now let us return to the question of attractive gravity. Observations indicate 
that gravity on the galactic scale is attractive (clustering, structure formation 
etc). By the existence of stable circular orbits we already know that the partic;les 
are being accelerated towards the galactic center. This can also be seen by 
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studying the geodesic equation for a test particle that has been "placed" at 
some radius ro, which yields 
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"2ri-' 
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dr 
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The quantity in square brackets is g"^^ and therefore must be positive. There- 
fore this expression is negative and it can be deduced that particles are attracted 
towards the center. Detailed discussions on this quantity may be found in Ford 
& Roman (1996) and Lobo (2008). A related issue is the following. In the New- 
tonian limit the gravitational energy of a spherically symmetric, gravitationally 
attractive system, with sufficient fall-off properties is negative. One could po- 
tentially argue that the Ad gravitational effects due to higher dimensional Weyl 
stresses on the brane should not be subject to the same criteria as in true four- 
dimensional gravity. However, from the point of view of the Newtonian limit 
on the brane, as well as CMB studies, it is desirable that the contribution 
from the Weyl stresses not be too large away from the galaxy. The question 
then becomes one of how to quantify such an effect in the relativistic case. 

We choose as our quantifier the following: The total gravitational energy Eg 
in the halo region must be negative (Misner, Thorne & Wheeler 1973; Lynden- 
Bell, Katz & Bicak 2007), where by gravitational energy we mean the quantity 
defined below. It is necessary to explicitly verify it because a positive energy 
density does not always lead to attractive gravity (Nandi et al. 2009). Hence 
we shall calculate the total gravitational energy Eq between two arbitrary but 
fixed radii with r > (we have excluded the origin since the solution is not 
valid at small r because of singularity there, and the effects of the core region 
are manifest in the constants in the solution) . This quantity is given by (Misner, 
Thorne & Wheeler 1973; Lynden-Bell, Katz & Bicak 2007): 



rr2 ^ 

Eg = M-Em = ^'Jt [1-A3]pr^ 

Jri 



dr 



(23) 



whore I—A2 < by definition (proper radial length is larger than the Euclidean 

length) and 

/•r2 

M = 4tt pr'^dr, (24) 

J ri 

Em is the sum of other forms of energy like rest energy, kinetic energy, internal 
energy etc. The dark radiation energy density (12) can be rewritten as 
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(25) 



For very small values of D and for < ^ ^ 2, it follows from the integration in 
Eq.(24) that Eq < for arbitrary r2 > ri > 0, indicating that gravity in the 
halo is indeed attractive. It is observed that the galactic halo extends as far as 
bOOKpc (Sahni 2004). Thus as an illustration we might take in units of Kpc, 
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ri = 100, / = 10~^, then it follows from Fig. 2 that Eq < for the upper limit 
r2 < 500. We recommend a value oi D < 0.00001, because it optimally yields 
reasonable behavior for all relevant quantities in the entire range of the halo 
6 OBSERVATIONAL CONSTRAINTS 

We wish to contrast the brane world model with the observational parame- 
ters in the galactic halo. For the ease of calculation, we shall rewrite the metric 
(5) in the form 

ds^ = -e2*('-) dt^ + + {de^ + sin^ Od^^ ) (26) 

r 

where, for the solution under consideration, the metric functions are specifically 
given by 

$(r) = i[lnBo+nnr] (27) 

m{r) = ^[l-A-\r)]. (28) 

These functions are not necessarily the same as the potential and mass functions 
obtained from observations. Note from equation (28) that A > I, which is a 
basic condition to be satisfied by a valid solution. In the first post Newtonian 
approximation, the gravitational potential $(r) is given in general relativity by 
the equation (Misner, Thorne & Wheeler 1973): 

» Rtt » Mp + Pr + 2pt) (29) 

which reduces to the equation for Newtonian i>jv where 

V24'Ar = 47r/? (30) 

only if the pressure contributions are negligible in comparison to energy density. 

The above general relativistic $ appears in the wavelength shifts z± of an 
emission line of a massive probe particle for an edge-on galaxy is (Nucamendi, 
Salgado & Sudarsky 2001; Lake 2004) 



1 

■z± 



^1 - r$'(r) 



(31) 



where b is the impact parameter and ± signs refer to approaching and receding 
particles. The above equation approximates to (Faber & Visser 2006) 

4 « r$'(r). (32) 

Therefore, the usual techniques for obtaining the potential for rotation curve 
measurements yield a "pseudo-potential" 

$RC = $ 7^ $jv. (33) 
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Prom Eq.(29), one then obtains a "pseudo-mass" as 



nine 



r^$'(r) « 47r j {p + Pr + 2pt)r^dr. (34) 



The right hand side expression comes from the first post Newtonian order in the 
weak field hmit and we call it MpM- We have so distinguished it because r^<l>'(r) 
comes entirely from pseudo potential $ whereas Mp^ is an integrated quantity 
that requires the knowledge of the equation of state. The pseudo mass torc 
reduces to the Newtonian mass M(r) [see Eq.(24)] when pressure contributions 
are negligible in comparison to energy density. 

Photons can also be regarded as probe particles as they sense the gravita- 
tional field in the halo during their travel to the observer. The effect of gravity 
on photon motion can be calculated in terms of a refractive index n(r). (See 
also Boonserm et al. 2005.) Por the Schwarzschild gravity, the exact geodesic 
equation, including the equation for Shapiro time delay, of massless particles 
was expressed in terms of n(r) by Nandi & Islam (1995) via the idea of optical- 
mechanical analogy. This idea has been further extended by Evans et al. (2001) 
to arbitrary spherically symmetric spacetimes and by Alsing (1998) to rotating 
spacetimes. The geodesic motion of both massless and massive particles could 
be exactly expressed in terms of a single generalized refractive index N = v^Vjc 
where V is the 3— velocity of the particle. For light motion, V =- c/n. 

Por unspecified metric functions <i>(r) and m{r), Faber and Visser (2006) 
argue that n(r) = 1 — 2$iens + in which they define the lensing pseudo- 

potential as 

$(r) 1 f m(r) , 
$iens = + TT / -^dr. (35) 



Another pseudo-mass miens obtained from lensing measurements has been de- 
fined as (Faber & Visser 2006) 

miens = ^r2$'(r) + ^m{r). (36) 
The first order approximations of Einstein's equations yield 

Pi^) i^[2<ns(0 - m'^c(r)] (37) 

inr^iPr + 2pt) « 2[m'^c{r) " <„s(0] (38) 

where the right hand sides denote pseudo-density and pseudo-pressures. Fur- 
thermore, Paber & Visser (2006) defined a dimensionless quantity 

, ,M - + ~ ^ m^^c - miens (^a) 

where the pseudo quantities on the right hand side of Eqs.(37)-(39) determine, 
respectively, the observed density, pressure and equation of state. Using the 
specific metric functions (27, 28) in Eqs.(35)-(39), we can calculate the first 
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to 



order quantities $rc, mRc, *iens, mens, 2(m^c ~ "^lens) I J^'^'^ -m°' 

lens RC 

be expected from combined observations. 

Note that pseudo quantities are the observable quantities. The general pro- 
cedure to determine the metric can be stated as follows: Once one is able to 
observationally determine the profiles of pseudo quantities, one can work back- 
wards to find the corresponding metric functions. This is a kind of reverse 
technique observational astrophysicists use. If the observed pseudo profiles fit 
(up to experimental error) with the analytic profiles of a priori given metric 
functions, one can say that the solution is physically substantiated. Otherwise, 
it has to be ruled out as non-viable. In this sense, the observed pseudo profiles 
play the role of constraints on the possible metric solutions. 

For the present solution, we obtain the following constraint equations on 
$(r), m(r), the pressure profile and the equation of state: 



$RC = ^(InBo + nnr) 
rriRc = r^$'(r) = ^ 



(40) 
(41) 



{lnBo + l\nr) + 
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2(mRc - mJenJ = I 



1(1"^ + 31 + 6) 
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(42) 
(43) 



2 m'^c - ^'lens _ 1 (^^ + -I- 6/ + 8)ri+T - D{1''^ + 3/^ + Gl + 4) 
3 D{1^ + 3P + 61 + 4) + (Z2 + QI + 8)ri+T 

(44) 



If observations yield the same profiles as on the right hand sides, we can say 
that they support the metric functions $(r) and m(r) of Eqs.(27, 28). Assuming 
for the moment that these specific functions are really observationally valid, we 
ask what comparative features of the observables do we expect to see. These 
can be visualized by first comparing the functions <I>rc with $iens- We plot 
their difference in Fig. 3 taking the numerical values of constants mentioned 
earlier and observe that despite their varying expressions, the actual difference 
$RC — $iens is negligible, of the order of only ~ 10~^. Thus rotation curve and 
lensing measurements lead to practically indistinguishable inferences about the 
pseudo potential, at least in the first order. Similarly, Fig. 4 shows variation 
with distance of Mpj^{r) and of different pseudo-masses torc and miens- All the 
masses share the common feature that they increase with distance r, though 



3 2m(^„^ - m'^c 

^ l{l + l) 
~ l + A ■ 
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they begin to slightly differ in the far field. To ascertain the role of pressures 
from the observational standpoint, one has to look at Fig. 5 where we find that 
the difference m'^Q — mj^^^g is of order ~ 10"'', to first order. To ascertain the 
equation of state from observational parameters, we come to Eq.(44), which is 
plotted in Fig. 6. One immediately finds that w — > 1/3. The closeness to 1/3 
actually depends on the value of D, the smaller its value the better the closeness. 
The exact equation of state following directly from the Eqs.(12)-(14) is 

Pr + 2pt = p (45) 

which implies from Eq.(39) that 

= i. (46) 

(Had the dark radiation been a perfect fluid, we could have p = ^p, the usual 
equation of state for radiation.) Noting that U = p and using Eq.(45) in Eq.(34) 
we find that the mass in the first post Newtonian approximation becomes 

Mpjv(r) =4n j{p + Pr + 2pt)r'^dr = 2M(r) (47) 

which is twice the Nciwtonian mass M(r) defined by Eq.(24). All the above 
indicate that the galactic halo modelled by the present solution shows somewhat, 
but not quite, Newtonian features (pressure contribution is not negligible as it 
is proportional to energy density). The profiles of (pseudo) potentials, masses 
and the anisotropic equation of state as expressed in Eqs.(40)-(44) above are the 
ones one would expect to find from the combined measurement of the pseudo 
quantities, if the solution (27, 28) is really a viable one. 

7 COMPARISON WITH A SCALAR FIELD MODEL 

Scalar fields are an integral part of reality though they yet lack any observa- 
tional evidence. We choose the well discussed scalar field model with potential 
V{r) proposed by Matos, Guzman & Nunez (2000) and compare it with the 
present brane world model. Using the flat rotation curve condition, they obtain 
the full solution as (we distinguish their quantities by tilde): 

B{r) = Bor' (48) 

^(^^ = 4 + i^(4-"/>-('+^) 

Hr) = J^l^r + (l)o (50) 

V OTT 

where D is an arbitrary constant of integration. For simplicity, they take D = 0, 
but this assumption makes A < 1 [see equation (28)]. Because of this, their 
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expressions for density and pressures pressures lead to different conclusions than 
those in the brane-world model considered here. For example, the expression 
for density exhibit p < 0, meaning violation of Weak Energy Condition (WEC) 
and furthermore it leads to 2 < —1, meaning repulsive gravity in the halo, 
contradicting observational facts. Therefore it is necessary to re-calculate the 
relevant quantities with D ^0. 

We find the density and pressure profiles in the rest frame of the fluid as 



We wish to emphasize here that the role of non-zero value of D is crucial not only 
for avoiding repulsive gravity but also for arriving at a correct conclusion about 
the relative strengths between pressure and density. For instance, let us take 
D = 1. In the distant halo region, we can take, say, r ^ 100 — 300 Kpc and with 
I ~ 10~^, we find the numerical values to be p ~ 10^^ and Pr ~ 10~®, which 
means that they are of the same order. But Pr+2pt 10~^^ => pr+2pt ~ 10~^p, 
which indicates that total pressure is roughly one hundred times less than the 
density. All these clearly go against the conclusion of the authors. However, if 
we take D = 0.00001, we find that Pr + 2pt 10^ jo. If we keep on decreasing the 
value of D further (but never exactly to zero for reasons stated above) we see 
that the total pressure dominates more and more over density so that the system 
becomes indeed non-Newtonian as claimed by Matos, Guzman & Nunez (2000). 
Such minuscule values of D are remarkably similar to those recommended for 
our solution (see Sec. 5). 

The next question is how far can we decrease D7 We notice the following 
interesting scenario: When D = 10~^, we find pr + 2pt = 9 x lO^p, which leads 
to LJ = ^''^j^* = 3 X 10^. This is the extreme non-Newtonian model possible 
in the scalar field model. On the other hand, if D = 10~^, we find that w > 
up to r = ro = 200 Kpc (attractive gravity) and becomes uj < —1 after r = tq 
(repulsive gravity). At r = ro, there is a singularity in w. When D = 10~^, 
we find that p < 0, to < —1, exhibiting the characteristics that follow from the 
same choice D = 0. Therefore we conclude that the limiting value of D is 10~^. 

The equation of state is anisotropic, as is evident from equations (54)-(55), 
a feature shared also by the brane world solution. It can be verified that p>0, 
p + Pr > 0, p + pr + 2pt > for D > 10~^, so we can say that the halo matter 
is not exotic because the WEC and NEC are satisfied everywhere. Therefore, 
we expect an attractive halo. To confirm it, again we follow the prescription in 
Lyndcn-Bell, Katz & Bicak (2007), and find that the total gravitational energy 
is indeed negative: 




(52) 



(53) 



(54) 




(55) 



13 



due to the fact that p > 0, 1 — < and r2 > ri. 

Certainly, the extreme scalar field model corresponding to D = 10"'' is 
highly non-Newtonian, that is, Pr + 2pt ~ l{fip. As a result, Eq.(30) leading 
to a purely Newtonian definition of mass M{r) as in Eq.(24) does not apply. 
However, incorporating the pressure contribution, the dynamical mass in the 
first post Newtonian order is 

Mpjv(r) =4tt J{p+pr + 2ptydr = 10''M(r), (56) 

which clearly reflects the non-Newtonian nature of the model in terms of masses. 
The observable pseudo quantities for this extreme case work out to 

mRc(r) = ^ « (57) 

„ , , l(l^ + l-4)r , , 

mie„s(r) « \(^i2_^) ^ 10"^ (58) 

2(mkc - miens) « « 10"' (59) 

~. ^ ^ 2 <c-<ns = l{l'-l-4)r'+^-D{l' + l'-Al-4) ^ 
~ 3 2m(^„^ - m'^c ^[^(ZS + /2 _ 4; _ 4) + ^2^2+!] 

(60) 

for our chosen range, r ^ 100—300 Kpc. Note that if we straightaway put D = 
in equation (60), we get u){r) < —1, conveying a completely different conclusion. 
Thus while the numerical values of observable potentials and masses behave like 
those in the brane solution, it is only the equation of state (61) that differs widely 
because a; has a value ^ 10^ compared to u = ^ [cf.Eq.(47)]. 
8 CONCLUSIONS 

In the foregoing, we first discussed the motivation for considering higher 

dimensions, which probably distinguishes the solution (27,28) as a more inter- 
esting search for the halo model. It is the 3— brane that we experience and 
therefore the bulk contributions U and P are translated into brane quantities 
via Einstein equations, which provide a way to have insights into the nature 
of the galactic fluid. The present model shows that the halo does not have a 
perfect fluid equation of state but NEC is preserved, meaning that the fluid is 
not exotic. It was shown that circular orbits in the halo are stable and that the 
dark radiation is attractive in nature, as it must be. Thus the solution satisfies 
two crucial physical requirements: Circular orbit stability and attractive gravity 
in the halo. The next task was to derive the constraints on the solution imposed 
by combined observations of rotation curve and lensing. 

The analyses by Nucamendi, Salgado & Sudarsky (2001), Lake (2004) and 
Fabcr & Visser (2006) have been carried out without reference to any specific 
form of metric functions. The specific functions (27,28) yielded the right hand 
sides of Eqs.(40)-(44) in terms of parameters to be measured by a combination of 
rotation curve and lensing measurements. These are effectively the constraint 
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equations on the chosen solution. If combined observations turn out to taUy 
with the behavior predicted by Eqs.(40)-(44), then the brane world solution 
(27, 28) can be said to be supported by observation. Until that happens, the 
solution would remain largely an academic curiosity. 

As discussed in the introduction, the increase of mass linearly with r comes 
from an elementary Newtonian argument to explain the observed fact of flat 
rotation curves. A key question still remains: How much of pressure contri- 
bution is there in the making of that mass, that is, is it just the Newtonian 
M(r) of Eq.(24) or Mpjv(r-) of Eq.(34)? The solution by Rahaman et al. (2008) 
does exhibit a linear increase of the Newtonian mass M{r). Fortunately, even 
if pressures are included, we find MpN{r) = 2M{r) [cf.Eq.(47)], that is, they 
are of the same order showing that the theoretical linear mass increase is pretty 
consistent with both the definitions within the present model. Such an in- 
crease is also supported by observable pseudo masses, as expected. However, 
this is just one aspect of the solution shared by many other models too. The 
distinguishing^ features lie elsewhere. For instance, we saw that in the scalar 
field model Mpjv(r) = lO^M(r) [cf.Eq.(56)]. As emphasized earlier, such dis- 
tinguishing features will have to be determined only by the detailed analyses 
of data profiles of potential, pressure and most importantly, the equation of 
state obtained through the combined rotation curve and lensing measurements. 
The bottom line is that the present solution depicts the halo as a mildly non- 
Newtonian model (pressure contribution is of the same order as energy density) 
in contrast to the highly non-Newtonian scalar field model. Although both the 
models share many properties including an anisotropic non-perfect fluid equa- 
tion of state, sharp differences appear in the computation of Mp^ir) and in the 
equation of state characterized by w(r). Observations to date do not seem yet 
conclusive enough (for a discussion, see Faber & Visser 2006) to tell us which 
one, if any, is closer to reality among all the competitive halo models. 
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